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We consider the non-equilibrium dynamics of topologically ordered systems driven across a contin- 
uous phase transition into proximate phases with no, or reduced, topological order. This dynamics 
exhibits scaling in the spirit of Kibble and Zurek but now without the presence of symmetry breaking 
and a local order parameter. The late stages of the process are seen to exhibit a slow, coarsening dy- 
namics for the string-net that underlies the physics of the topological phase, a potentially interesting 
signature of topological order. We also note a time dependent amplification of the energy splitting 
between topologically degenerate states on closed manifolds. We illustrate these phenomena in 
the context of particular phase transitions out of the abelian Z2 topologically ordered phase of the 
toric code/Z2 gauge theory, and the non-abelian SU(2)fc ordered phases of the relevant Levin- Wen 
models. 



I. INTRODUCTION 

This paper lies at the intersection of two interesting 
streams of contemporary research: the study of the non- 
equilibrium dynamics of quantum systems, and the study 
of topologically ordered phases of matter. Consider driv- 
ing a topologically ordered system through a phase tran- 
sition to a topologically trivial, or relatively trivial, phase 
by changing some parameter in the Hamiltonian at a 
slow, but finite rate. This work investigates the ensu- 
ing non-equilibrium dynamics resulting from this "trans- 
critical protocol"^, with particular emphasis on univer- 
sality and the non-Landau character of the transition. 

The Kibble-Zurek (KZ) mechanism^] 

is a scaling the- 
ory of the defects generated by slowly cooling a classical 
system through a continuous symmetry-breaking phase 
transition. Kibble originally formulated this problem 
as a cosmological theory of the phase transitions in an 
expanding universe, and Zurek later applied it to con- 
densed matter systems. The central physical insight is 
that sufficiently distant parts of the system settle into in- 
dependently seeded local broken symmetry "directions" 
whence the mismatch must be accommodated by a finite 
density of defects. The mechanism has since been gen- 
eralized to quantum phase transitions^^ and has been 
experimentally investigated in a number of systema^^, 
although decisive confirmation of the predicted critical 
scaling of defect densities is still lacking^. Polkovnikov 
and co-workers have also studied the scaling theory of the 
excess heat density an d the interplay between the ramp 
velocity and finite siz d 12 * 13 l For a recent review of the 
broader context of these developments in the study of 
non-equilibrium q uantum phenomena, see Ref. [14] . 

Few obser vable d 1 3 1 1 5 1 1 £l other than the universal non- 
equilibrium defect density have been studied in the 
KZ problem, now defined more broadly as the non- 
equilibrium temporal evolution of a system in the vicinity 
of a critical point. Recently, Chandran et. alP^have sys- 
tematized the universal content in the KZ problem in a 
scaling limit, and written non-equilibrium scaling func- 



tions for all physical observables in this limit. The scaling 
functions describe the entire time history, and asymptote 
to equilibrium and coarsening scaling regimes in the ap- 
propriate limits. The universal content depends only on 
the pairing of the dynamical universality class of the crit- 
ical point, and the particular protocol through parameter 
space. 

The KZ scaling theory also provides an elegant frame- 
work within which to investigate ramp dynamics of phase 
transitions beyond those of traditional symmetry break- 
ing. These transitions could involve the destruction of 
the topological order of states of matter like spin liq- 
uids and the fractional quantum Hall phases. Topo- 
logically ordered phases are not locally distinguished by 
any order parameter, but are characterized by emergent 
gauge fields and fractionalized excitations. Their non- 
local structure makes them particularly robust to local 
perturbations and well-suited to perform "topological" 
quantum computation^. Despite their robustness, a 
strong enough perturbation can drive a transition from 
a topologically ordered phase to a trivial or relatively 
trivial (i.e. one with a smaller gauge group) phase. For 
appropriately selected perturbations this transition will 
be continuous. 

This paper addresses the KZ problem when a system 
is driven from a topologically ordered phase to a prox- 
imate trivial or relatively trivial phase. We do so for 
a class of topological phases that possess lattice realiza- 
tions where the gauge degrees of freedom are mani fest: 
these are the toric code/lattice Z 2 gauge theorySED and 
the string-net models of Levin and WerPSl that realize 
doubled non- Abelian Chern-Simons theories. By a com- 
bination of duality and perturbative arguments, we show 
that KZ scaling in the generalized sense of Ref. [T] holds 
for various observables even though the canonical KZ sig- 
nature of a density of topological defects is not meaning- 
ful. We further provide strong arguments that the late 
time dynamics in the scaling regime exhibits a slow coars- 
ening of the string-net that is condensed in the starting 
topologically ordered state. To our knowledge, this is the 
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first treatment of a quantum coarsening regime in the dy- 
namics of an isolated, quantum system. As the extended 
string-nets are central to the topological character of the 
starting phase^, their slow decay outside the phase is a 
(potential) signature of the physics of the parent phase. 
In the same spirit, we find that the exponentially small 
energy splitting of topologically degenerate ground states 
on closed manifolds can be amplified in the KZ process 
to a splitting linear in the system dimension. The re- 
striction to the scaling limit always brings simplification 
as particular gapped degrees of freedom are, at worst, 
dangerously irrelevant. That is, they do not affect the 
scaling regime but do alter the asymptotically long time 
behavior of the KZ process. 

The most relevant precursor to our work is found in 
the cosmo logy l iterature in the papers of Rajantie and 
Hindmarsr i 22 l 23 [ They studied the non-equilibrium dy- 
namics of ramps through the finite temperature phase 
transition in the non-compact Abelian Higgs model; their 
protocol moves between the gapless Coulomb phase with 
gapped matter to the fully gapped Higgs phase. However, 
there are three important differences. First, their work 
involves finite temperature in an essential way. The zero 
temperature limit of their protocol would involve exciting 
the system in the gapless phase even before the transi- 
tion is reached. Second, the non-compactness of their 
gauge field makes the physics of their Higgs phase quali- 
tatively different from that of the compact gauge models 
considered by us. This difference is quite visible in our 
choice of observables. Third, we work on the lattice in 
the "electric flux" representation, a natural choice in the 
condensed matter setting, while they work in the contin- 
uum with the vector potential. Thus our discussion of 
string-net coarsening has no analog in their formulation. 
In the condensed matter literature, ramps across topo- 
logical transitions^ in (1+1)D and sudden quenches in 
one of our model sy stems , the perturbed toric code, have 
been studied before^^U. The sudden quenches are in a 
completely different limit from the slow ramps we study 
here as they inject a large amount of energy into the sys- 
tem. Finally, low temperature spin ice exhibits topologi- 
cal order in a classical limitP^l and its dynamics following 
quenches is dominated by monopoles of the gauge field. 
However in this case, non-universal lattice effects turn 
out to dominate the long time behavioral . 

We turn now to the contents of the paper. We begin in 
Sec. [H] by reviewing the phase diagram of the Z 2 gauge 
theory coupled to matter and describing the transitions 
out of its topologically ordered phase. Readers literate 
in the canon of topological phases can skim this sec- 
tion for our notation. Section III describes the KZ ramp 



across the pure matter sector of the Z 2 theory which has 
a conventional symmetry-breaking transition. This sec- 
tion contains a new analysis of coarsening in the (2+l)D 
transverse field Ising model and summarizes our previ- 
ous understanding of KZ. In |IV| we discuss our results 
on the KZ scaling functions and string-net coarsening for 
ramps in the pure gauge sector of the theory, which has a 



confinement transition without a local order parameter. 
We then generalize the scaling theory to a ramp across 
an arbitrary point on the critical line in the phase dia- 
gram in [Vj In |VI[ we discuss the energetic amplification 
of the topological degeneracy in the deconfincd phase by 
the KZ protocol. We then turn to generalizations of these 
results to phases with non- Abelian topological order in 
Sec. |VII| specifically, we discuss a particular transition 
from the SU(2)k ordered phases. We end with a discus- 
sion of generalizations to other theories. 



II. REVIEW OF THE PHASE DIAGRAM OF 
THE Z 2 GAUGE THEORY 

The phase diagram of the (d+1) dimensional Z 2 gauge 
theory with matter^ contains a topologically non-trivial 
(deconfined) phase and a topologically trivial (confined- 
Higgs) phase. The topological order in the deconfined 
phase is described by the BF theory^. We work in d = 2 
for which the Z 2 theory is precisely (the topologically 
ordered) Kitaev's toric code with perturbations^. We 
start by reviewing the key features and the excitation 
spectrum of the toric code. We then discuss two pertur- 
bations that drive a continuous transition to a topologi- 
cally trivial phase. 



A. The perturbed toric code 



The toric cod c! 19 * 31 J is defined in terms of spin- 1/2 de- 
grees of freedom that live on the links I of a 2D square 
lattice: 

H TC = —K ^ Bp — r m ^ A s 
p s 

^En^- r ^E n «? w 

p ledP s i.sedi 

where the A s and B P are "star" and "plaquette" opera- 
tors, s and P denote the sites and elementary plaquettes 
of the lattice, while dP and dl are the boundaries of 
plaquettes and links. H TC can be rewritten as a gauge 
theory with matter by identifying the tr; variables as the 
gauge degrees of freedom, and introducing new spin 1/2 
'matter' variables, r s , on the sites of the lattice. Upon 
restricting the expanded Hilbert space to the 'physical' 
subspace of gauge-invariant states 



(2) 



i.sedi 



the toric code Hamiltonian ([lj is equivalent to the gauge- 
invariant Hamiltonian: 



P iGdP s 



(3) 



Note that G s defines a set of local symmetries at each 
site since [H , G s ] — V s. 
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FIG. 1. A section of the toric-code lattice with operators in the magnetic (left) and electric (right) bases. Matter (r) and gauge 
(a) variables are located on the sites and links respectively. The W m operator |4| flips a string of a z variables and creates a 
pair of magnetic vortices at its endpoints, while the W e operator flips spins in the x basis to create a pair of electric charges 
linked by electric flux. 



In the x basis of the spin operators, it is useful to think 
of r and a as the electric 'charges' and 'fluxes' in the 
theory respectively: t£ = —1 (+1) if an electric charge 
is present (absent) at site s, while af = —1 denotes the 
presence of electric flux on link I. In this language, we 
recognize the gauge-invariant condition ([2]) as the lat- 
tice Z2 version of Gauss's law. In the conjugate z basis, 
the operator Bp = JlzedP °f measures the magnetic flux 
through the plaquette P. 

The model is exactly solvable as both terms in H com- 
mute with each other and the gauge constraint G s . The 
ground state is charge-free and vortex-free: a simulta- 
neous +1 eigenstate of t£ and Bp for all s,P. As the 
ground state is free of charge, it is a loop gas of elec- 
tric flux. That is, it is an equal amplitude superposi- 
tion of configurations where links with af = —1 form 
closed loops. The degeneracy of the ground state mani- 
fold depends on the topology of the lattice; on the torus, 
it is four-fold degenerate. The four ground states can- 
not be distinguished locally. They are labelled by the 
eigenvalues ±1 of the non-local Wilson loop operators 
W nc = Y\i eC of along the two distinct non-contractible 
loops C nc on the direct lattice. 

The elementary excitations of the model are gapped 
and are of two types: e and m. e denotes the presence 
of electric charge on site s, while m is a magnetic vortex 
on plaquette P characterized by Bp = — 1. e and m are 
individually bosonic, but have mutual semionic statistics. 
The non-local string operators 

W e (s >S ')= J] r s Vfr*, W m (s,s')= J] of 
s.s'edC s.s'edc 

(4) 

defined respectively on the curves C and C on the direct 
and the dual lattice, create a pair of electric charges and 
vortices at their ends as shown in Fig. [l] 

H is robust to small local perturbations and ex- 
tends to a topological phase. Nevertheless, a strong 



enough perturbation will eventually drive a transition 
into a trivial phase. The toric code Hamiltonian per- 
turbed by transverse fields is H = Htc — J2i + J°i 5 
both perturbations drive continuous transitions to trivial 
(spin-polarized) phases when made large. In the gauge- 
invariant formulation of the Z 2 gauge theory with matter, 
the perturbed Hamiltonian takes the form 

p s 1 s edi 1 

(5) 

The phase diagram of this theory was explained in de- 
tail in the seminal paper by Fradkin and Shenker^Sl and 
has more recently been confirmed in several numerical 
studies^HUl. We will now briefly review this model in 
different parameter regimes. 

B. Pure matter theory (T = 0) 

For r = 0, the gauge degrees of freedom are static 
and frozen into a vortex-free configuration in the ground 
state sector. It is therefore convenient to diagonalize H 
in the gauge- variant subspace where erf = I for all I and 
project the eigenstates to the gauge-invariant subspace 
afterwards. The Hamiltonian then maps to the (2+l)D 
transverse field Ising model (TFIM) for the matter spins: 

H TPIM =~JJ2 T s' T s r M E T s- ( fi ) 
< SS '> s 

On tuning J, the TFIM undergoes a conventional 
'Higgs' phase transition from a paramagnetic phase to 
a symmetry-broken ferromagnetic phase. In a comple- 
mentary view, the static electric excitations e defined at 
the toric code point (T = J = 0) acquire dynamics when 
J =/= and eventually condense at a critical value of J. 
The transition is in the 3D Ising universality class and is 
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detected by the local order parameter, (r z ) in the gauge- 
variant subspace. 

After projection, the state with {rj ; z = 1} is the vortex- 
free configuration in the topological sector defined by the 
Wilson loop W nc = 1 for both non-contractible loops 
on the torus. This choice maps to a TFIM with peri- 
odic boundary conditions in both directions of the torus; 
the remaining three topologically inequivalent vortex- 
free configurations generate TFIMs with different bound- 
ary conditions (periodic-antiperiodic etc.). The four- 
fold degeneracy of the topological phase vanishes in the 
Higgs/ferromagnetic phase. 



C. Pure Gauge Theory (J = 0) 

In this case, matter is static. The ground state is in 
the charge-free sector (t£ = 1), and the Hamiltonian for 
the gauge variables in this sector is: 



II- 



(7) 



When T/K is small, the gauge variables are weakly fluc- 
tuating and the elementary excitations are well-described 
as vortex pairs. At some critical T/K, the vortices con- 
dense and the gauge variables strongly fluctuate past this 
point. This transition cannot be diagnosed by a local or- 
der parameter. Instead, the vortex condensate phase is 
marked by the vortex pair creation operator, (W m ) =/= 
for vortices separated by long distances. In the conju- 
gate electric field basis, flux loops become costly as T is 
increased; hence the transition is from a topological loop 
gas phase at the toric code point to a phase in which flux 
loops become confined. 

In a different language, the transition is understood 
as a deconfinement-confinement transition for the static 
electric charge^ and is diagnosed by the free energy cost 
of creating a pair of (infinitely separated) charges. The 
cost is finite in the deconfined phase, but infinite in the 
confined phase, and is equivalent to the change in behav- 
ior of the expectation of the contractible Wilson loop: 
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FIG. 2. T — phase diagram of the Z2 theory in d = 2 
dimensions. The matter and gauge axes are dual, and the 
Higgs and charge confined phases are smoothly connected. 



D. The full phase diagram 

The full T = phase diagram of the Ising gauge theory 
in (2+l)D is shown in Fig.pl Fradkin and Shenker^Slhave 
shown that the confinement/Higgs transitions are stable 
on moving away from the pure gauge/matter axes. Fur- 
ther, the Higgs and confined phases are smoothly con- 
nected. However the diagnostics previously discussed, 
like the Wilson loop, no longer differentiate between the 
two phases. In a recent paper^, Gregor et. al. have 
shown that an appropriate ly def ined line tension, related 
to the Fredenhagen MarcrP^' order parameter studied 
by lattice gauge theorists, can be used to diagnose the 
transition everywhere in the phase diagram. We will use 
this quantity in combination with the topological ground 
state degeneracy to study ramps across generic points on 
the critical line in the phase diagram. 



W(L) 



wee 



(8) 



from a perimeter law (W(L) ~ exp(— L)) to an area-law 
(W(L) ~ exp(— L 2 )). C is a contractible loop and L is 
its perimeter. 

In d = 2, the Z2 Ising gauge theory is self dua l 18 * 36 ! 
Thus the pure gauge theory also maps to a (2+l)D TFIM 
and the confinement-deconfinement transition belongs to 
the 3D Ising universality class. The details of this duality 
are explained in Ref. 1181 and have been summarized 
in Fig. [3] We emphasize that despite the duality, the 
transition is not described by a local order parameter, as 
will be even clearer in Sec. IVIII 



III. KIBBLE ZUREK I - RAMP ACROSS THE 
HIGGS TRANSITION 



We begin by reviewing the Kibble-Zurek (KZ) formal- 
ism for linear ramps (a Trans-Critical-Protocol in the 
parlance of Ref. H} across the conventional 3D Ising tran- 
sition along the pure matter line at T = in Sec. |III A| 
The late time evolution of the system is naturally de- 
scribed as "coarsening" : a dynamical process previously 
discussed only in classical systems. We take the first 
steps to apply these ideas to an isolated quantum system 
in Sec. ITTTF31 
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A. Scaling theory: Review 

The system is in equilibrium in the paramagnetic phase 
at t = — oo and is driven to the ferromagnetic phase 
by changing the transverse field Tm^) linearly: 8{t) = 
CTm(*) — ^Mc)/Tmc — t/ T - The critical point (CP) is at 
6 = and r is the ramp time. The response of the sys- 
tem to slow ramps is characterized by three chronological 
regimes: adiabatic evolution at early times, 'critical' or 
diabatic evolution near the CP, and a late-time regime 
that we argue to be domain-growth. 

At early times, the system is far from the CP and 
evolves adiabatically. Critical slowing down implies that 
the instantaneous correlation time diverges as £*, ~ — 
Ta/cI _lyz near the CP and adiabaticity must break down 
before the CP is reached, v and z are respectively the 
correlation length and dynamic exponent with v = 0.627 
and z — 1 in this case^l . The system falls out of equilib- 
rium at the KZ time, t = —t K , when the time remaining 
to reach the critical point, t K , becomes equal to £ t : 



U 



(9) 



At t = —t K the evolution of the system become diabatic 
and, to zeroth order, the system remains frozen until it 

emerges on the other side of the CP at t = t K . The 

l/z 

KZ time, t K , defines a KZ length, l K = which is 
the correlation length at the time the system falls out of 
equilibrium. Finally, for t 3> t K , we get coarsening. 

Recently, Chandran et. alP formulated a scaling limit 
in which the physics described above becomes universal 
for a given pairing of a critical point and a ramp proto- 
col. This limit is defined as r — > oo with time and length 
scales measured in units of the diverging scales t K and 
l K . As 5{t K ) — > in this limit, the out-of-equilibrium 
response of the system is completely controlled by the 
critical point. The content of the scaling theory is not 
just the critical exponents, previously discussed by Kib- 
ble and Zurek, but also scaling functions for various phys- 
ical observables. For example, the equal-time two-point 
correlation function of the order parameter defined as: 

(r:(t)r:,(t)) T ^G TT (\s-s'\,t;T), (10) 

has the KZ scaling form 



lim l 2 K A G TT (x,t) = G TT 



1 ' t 



(11) 



where t — >* oo is defined as the limit r 



with 



and held fixed, x is the distance between sites s and 
s' , and A is the scaling dimension of the operator t z . A 
is equivalent to /3/v, and is numerically found to be 0.518 
in the (2+l)D TFEVP . 

The scaling function should asymptote to the correct 
equilibrium form in the limit t/t K — oo with x/£(t;r) 
fixed 



where x and t are defined as x/l K and t/t K respectively, 
Q eq is the equilibrium scaling function, known to decay 
exponentially in the Ising model, and xt v — x/£(t). 

A full description of the diabatic regime, previously 
termed 'critical coarsening^, is much harder. The dif- 
ficulty surpasses static computations in the analogous 
quantum critical regime as it involves real time. 

As the scaling content at late times has not appeared 
in the literature before, we devote the next subsection 
to it. We can also investigate the scaling functions for 
thermodynamic quantities like the excess energy density 
above the ground state^, q, and the entropy densitjEHH 



q(*,T)~l-*t- 1 Q(t) 
S (t;r)~l- d S(t) 



(13) 
(14) 



Both quantities tend to zero as t/t K — > — oo when the 
evolution is adiabatic and become non-zero when the sys- 
tem falls out equilibrium. The evolution at late times is 
therefore best understood through the finite excess en- 
ergy density (q > 0) or finite temperature phase diagram. 



B. Coarsening 

We now address the late-time dynamics of the KZ 
ramp in three steps. First, consider a KZ ramp from 
the disordered to the ordered phase across the 2D Ising 
phase transition at finite temperature. This is almost 
the classic KZ mechanism. The lack of long range or- 
der in the system when it emerges into the ordered 
phase drives the subsequent process of domain growth 
called coarsening. At late times, the system is locally 
in equilibrium with long domain walls on the growing 
length scale, l co {t) » l K . The physics of this process 
is captured in a classical hydrodynamical theory on the 
coarse-graining scales l K , t K belonging to Model C in the 
Halperin-Hohenberg classification 4 -. In addition to the 
non-conserved order parameter field in the oft-studied 
Model A, Model C has a conserved energy density im- 
posed by the underlying quantum dynamics. Coarsening 
is then an entropic effect in which the decreasing energy 
density in the dilute domain wall network is compen- 
sated by the increasing bulk energy density. This be- 
havior has been numerically verified in a few microscopic 
modelsSS] 

Now consider the KZ ramp across the quantum critical 
point in which we stop the ramp at some t/t K = t s 3> 
1, while continuing to measure time and length on the 
scales set by t K ,l K . The superscript s denotes stopping. 
Unlike in the previous case where q (Eq. 



131 is finite 



g„ (x,i)~P^g:«(xi») 



(12) 



and set by the critical temperature, q here vanishes in 
the limit of an infinitely slow ramp. The justification 
for the emergent hydrodynamical theory is then slightly 
different. First, we show that the evolution for t/t K > t s 
is indeed in the finite-temperature ordered phase. To 
wit, note that the dominant contribution to q{t s ;r) is 
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from the defect density on the scale l K frozen in at t « 
t K . Assuming that these defects evolve adiabatically for 
t > t K , we may conservatively estimate g(t s ;r) to scale 
as the single-particle gap l/£ z . As promised, this density 
is much smaller than the instantaneous critical density, 
Qc ~ l/£, d+z , above which the system is disordered in 
equilibrium: 



q ~ e ~ (! I J 



(15) 



Once we stop in the ordered phase, the time-scale for 
local equilibration processes on the scale l K is set by t K 
and t s . Coarse-graining on these scales, we may once 
more employ Model C to describe the slower dynamics 
on the longer length scale of the frozen string defects^. 
Specifically, we then predict that the system obeys the 
dynamic scaling hypothesis, that is, that it looks self- 
similar on the scale of a growing length l s co (t\ r) as t — > oo: 



gi (x,i) ~ (t**) a " A 0£ {xi K /v co ), 

l/z d 



(16) 



where l s co {t; r) = l K 



and Zd = 2. 



The value of the dynam ic exp onent, Zd, quoted above is 
only known numerically^'^', is a scaling function 
that can be computed within Model C. 

Finally, in the original KZ problem (where we don't 
stop the ramp), the late time evolution is also in the 
finite-temperature ordered phase as the relation q/q c <C 1 
holds for every t/t K 3> 1. However, the continuously 
changing parameter in the Hamiltonian affects the late- 
time limit in two important ways. First, the characteris- 
tic size of the domains grows at a slower rate: 



^Co(^7 7") — lh 



t 



where 9 



Zd 



1 

Zd ' 



(17) 



In the (2+l)D TFIM, 6 = (1- v)/2 and is smaller than 
l/zd- This slowing down can only help in the justifica- 
tion of a hydrodynamical description. The second effect 
is that the single particle gap A grows as l/t K {t/t K ) uz 
at late times. This, however, increases various scattering 
times in the problem and the applicability of hydrody- 
namics here becomes a delicate affair. In Appendix [Aj 
we argue that the process that drives coarsening and in- 
creases entropy involves the interaction of the long do- 
main walls with the bulk quasiparticles within each do- 
main. As the bulk quasi-particles scatter off the walls 
parametrically many times before the system parame- 
ters are changed, coarsening can at least self-consistently 
be justified. We therefore conjecture that the two point 
function as t — > oo holding x/l co (t) fixed obeys dynamic 
scaling: 

g TT (x,i) ~P» A gz(xi K /i co ). (is) 

In this process, the entropy density increases weakly 
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FIG. 3. Table summarizing the duality between the pure 
Z 2 gauge theory (7} and the TFIM in d = 2. Dark and light 
lines denote the direct and dual lattice respectively. 



in time. The late time asymptotes reflect this: 



Q{t) ~ qo i vz + qi t" 

s{t) ~ So - Sl i~ e . 



(19) 



The leading terms in S would be present even if the evolu- 
tion were adiabatic. The sub-leading term is the thermo- 
dynamic signature of coarsening. From this point, every 
time we invoke results from coarsening, the reader should 
keep in mind the subtleties presented in this section. 



IV. KIBBLE ZUREK II - RAMP ACROSS THE 
CONFINEMENT TRANSITION 

We now ramp across the pure gauge theory Eq. Q by 
tuning r. In this case, the transition is from a topolog- 
ically ordered deconfined phase to a confined one, and 
there is no symmetry-breaking description in terms of a 
local order parameter. Nevertheless, we will now show 
that the KZ mechanism for Landau transitions discussed 
in the previous section can be generalized to these transi- 
tions. Additionally, the loops and strings characterizing 
the topological phase (string-nets) will coarsen. 

Our main tool is the duality in (2 + 1) between the 
pure gauge theory and the TFIM summarized in Fig. [3] 
Importantly for us, the presence of electric flux on a link 
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(of the direct lattice) maps to a domain wall between the 
TFIM spins (on the dual lattice), while the vortex opera- 
tor Bp maps to the dual transverse field. The duality also 
ensures that a finite temperature confined phase exists, 
and that coarsening is described by the hydrodynamics 
of Model C. 

For specificity, we begin the ramp at the deconfined 
toric code point in one of the ground state sectors. The 
ground state is a loop gas of the electric flux lines in the 
a x basis. By duality, these are the domain walls of the 
paramagnetic phase of the TFIM. The system falls out 
of equilibrium in the deconfined phase before it is taken 
through the transition, with a network of loops of mini- 
mum size l K . In the confined phase, flux loops map to the 
costly domain walls of the dual ferromagnetic phase. Post 
the diabatic regime in the confined phase, this network 
of loops (string-nets) is diluted (average size increases 
as l co ) as the system coarsens. More generally, we can 
imagine string-nets being diluted in a generic topological 
theory and we will show some examples of this in Section 

As before, energy conservation requires the decreasing 
electric flux density to be compensated for by an increas- 
ing bulk energy density. Essentially, the system arrives in 
the confined phase (which is a vortex condensate) with 
a greater electric field density and a smaller magnetic 
vortex density as compared to the instantaneous ground 
state. The subsequent evolution through coarsening in- 
creases the typical size of the electric flux loops to l C o(t), 
thereby decreasing the electric field density and increas- 
ing the bulk energy density of the vortex condensate. 

Next, the two-point correlator that detects long-range 
order in the dual TFIM, (f§?§,), maps to the vortex pair 
creation operator Q, (W m (s,s')) that detects vortex 
condensation. As the condensed phase is also a confin- 
ing phase for charge, a non-zero value of (W m ) for long 
strings detects charge confinement. The scaling form for 
(W m ) is given by Eq. (11), and its asymptotic behavior 



Its asymptotic behavior is: 



is identical to that of the two-point function discussed 
in the pure-matter theory. In particular, in the coarsen- 
ing regime, the dual TFIM is ordered on length scales 
less than l co (t). Correspondingly, (W m ) is also non-zero 
on scales shorter than l co (t) but decays exponentially on 
longer length scales. Thus, the non-local string operator 
(W m ) probes the crossover scale from confinement to de- 
confinement as a function of time. Fig. |4](a) shows the 
scaling for (W m ). 

Finally, we can consider an interesting observable that 
we did not discuss in the TFIM. This is the Wilson 
loop @, W(R,t;r) on a curve of radius R. Were the evo- 
lution to be adiabatic, W(R, t; r) for large R would obey 
a perimeter law when t < and an area law when t > 0. 
In the KZ scaling limit, the scaling of the Wilson loop 
takes the form W(R,t;r) ~ W(R,i), where R/l K = R. 



W(R,i) 



exp(-i?t v ), if<<C-l 
exp(-(i?^) 2 ), if t > 1 and R < i e 
exp(-(i?/F)), if t > 1 and R > t 6 



(20) 

These scaling forms follow simply from the picture of adi- 
abatic evolution when t <^ —1 and a growing length l co {t) 
separating confinement from deconfinement when t 3> 1. 
The Wilson loop therefore also probes the crossover scale 
from confinement to deconfinement as a function of time. 
Fig. Qb) shows the scaling of the Wilson loop. 



V. KIBBLE ZUREK III - RAMP ACROSS A 
GENERIC TRANSITION IN THE Z 2 THEORY 

We will now see how the discussions of the previous 
two sections can be generalized to ramps crossing any 
critical point in the full Z2 phase diagram. First, con- 
sider moving off the pure gauge line by introducing a 
small, but non-zero J. The coupling to matter is irrel- 
evant to the T = transition; hence, the confinement- 
deconfinement transition in Fig. [2] persists for non-zero 
J and remains in the same universality class. Since the 
gap to charge excitations does not close on making J 
non-zero, we can re- write the Hamiltonian as one with no 
gauge-matter coupling (to any fixed order in J) through 
a canonical unitary transformation. The transformation 
defines "dressed" charge and gauge operators — in the 
dressed variables, the ground state is charge-free and A c 
is the non-zero gap to charge excitations. 

As we heat the system in the process of the ramp, 
we also need to consider the finite temperature phase 
diagram when J / 0. Although the ground state is 
(dressed) charge-free, the excited states have an expo- 



nentially small density of dressed charge, e 



-A c /T e// 



at 



any effective temperature r e // corresponding to the ex- 
cess energy density q in Eq. ( 19 ). The presence of charge 



at finite temperatures is extremely significant for the late- 
time coarsening picture for two reasons. First, it destroys 
the finite temperature confined phase at any non-zero J, 
without which a coarsening description is not meaning- 
ful. Synergistically, a finite density of charge implies that 
the electric field lines naturally end somewhere. Thus, 
the pictures of flux- loops/domain walls coarsening are 
no longer sensible at the longest length scales. 

Fortunately, the Kibble-Zurek scaling limit saves us 
from the problems raised above. This is because the ratio 
q/q c goes to zero as t/t K — > 00 (Eq. (15)) or equivalcntly, 



the effective "temperature" , T e f f , computed from q goes 
to zero in the scaling limit. A c , on the other hand, re- 
mains finite. This implies that the ratio of the average 
distance between charges to the KZ length, e A "^ T " ff /l K , 
is formally infinite in the scaling limit. Thus, while the 
dressed charges modify the true long-time behavior by 
ending coarsening, the scales on which they do so lie out- 
side the KZ scaling regime: in this way, the coupling 




FIG. 4. (a) The scaling function of (W m ) as a function of the scaled vortex separation |s — s'\ at different scaled times. This 
operator is dual to the two-point correlator (t|t|/) in the TFIM. (b) The logarithm of the scaling function of the Wilson loop 
as a function of the scaled radius illustrating Eq. (201. The scales on the x-axis are £,(t;T)/l K = i ~ v and l C o{t)/li< = i° ■ 



to matter is a dangerously irrelevant variable in the KZ 
problem (in the scaling limit). 

While we can write scaling functions for dressed ob- 
servables, the results are not very elegant since the 
Hamiltonian dependent dressed operators are different 
at different points in time. A crisper solution is to 
use the li ne te nsiorP^/ Fredenhagen-Marcu (FM) order 
parameteJ 38 * 39 * alluded to previously. This is defined as 



Sec. II B the couplings were all ferromagnetic as of = 1 



R(L) 



W 1/2 (L) 
VW(Lj 



JSC a i 



(21) 



where C is a square loop of side L and C\i<i is the open 
rectangle of sides L and L/2 obtained by cutting C in 
half, s,s' are the endpoints of C\/2 and W is the con- 
tractible Wilson loop. 

As L — > oo, R{L) is zero in the deconfined phase and 
non-zero otherwise. In this way, R{L) acts as a test 
of long range "order" , and appropriately generalizes the 
two-point spin correlator G TT from the pure matter the- 
ory (10), and the vortex pair creation operator W m Q 



from the pure gauge theory. In fact, in the gauge-variant 
subspace {of = 1} on the pure matter line, R(L) exactly 
reduces to G TT . The scaling form and asymptotes of R(L) 
therefore follows from Eq. (11 1 and the discussion below 



it. Of course, by duality, an identical analysis can be car- 
ried out by perturbing away from the pure-matter line as 
long as we interchange the gauge and matter degrees of 
freedom. 



VI. AMPLIFICATION OF TOPOLOGICAL 
SPLITTING 



Coarsening in the confined phase is one signature of 
the topological nature of the deconfined phase. We now 
present a second signature along the pure gauge line. 

Recall that the nearest neighbor coupling on each link 
in the TFIM on the pure matter line is set by — J erf- In 



for all /. The projected eigenstates of the TFIM were 
therefore in the topological sector labelled by the eigen- 
values (+1,+1) of the non-local Wilson loop operators 
W nc along the non-contractible loops in the x and the 
y directions. The other three sectors contain a line of 
anti- ferromagnetic bonds and hence an extra domain wall 
along the x, the y, and both directions on the torus re- 
spectively. The four sectors are degenerate in the para- 
magnetic phase. In the ferromagnetic phase, the domain 
wall cost per unit length is given by the single-particle 
gap, A; we therefore expect the splitting between the 
sectors to scale as the gap times L. This splitting has 
been confirmed numericalljPS 

Now consider the KZ ramp along the pure gauge line, 
but starting with general superpositions in the ground 
state manifold in the deconfined phase. By duality, we 
conclude that there can be a time-dependent contribution 
proportional to A(t; t)L to the total energy in excess of 
the ground state in the confined phase at late times. This 
is in addition to the leading contribution independent of 
the initial state in the deconfined phase, q(t;T)L 2 , that 
we introduced in Eq. ( 13 ). Substituting the known time- 



dependences of these quantities at late times t/t K 1, 
the excess energy of the ramp scales as: 



E(t; T ) 



1 / L \ 2 / t x vz 



L ( t 



Ik U 



where c is a geometry-dependent constant set by the ini- 
tial state. The apparent mismatch of units is because 
the lattice spacing has been set to 1. Thus, an energy 
difference between the topological sectors of magnitude 
exp(— LA(t; t)) at early times t/t K <C — 1, can be ampli- 
fied to O(L) at late times. More strongly, when c ^ 
(true for almost all superpositions in the ground state 
manifold in the deconfined phase), a contribution at O(L) 
with the time-dependence above implies that the ground 
state of the phase at early times contains extended string- 
like degrees of freedom, that is, it is topologically ordered. 
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VII. EXTENSION TO GENERALIZED 
LEVIN- WEN MODELS 

In this section, we will generalize the Kibble-Zurek 
problem to transitions out of a broader, non-Abelian, 
class of topological phases in lattice spin models of the 
Levin- WerP^ type. We will restrict our discussion to the 
subset of SU(2)fc models whose topological order is that 
of a doubled, achiral, Chern-Simons theory with gauge 
group SU(2) and a coupling constant of ft in appropriate 
units, though the construction of Ref. |2Q]is more general. 
We will also restrict to particular transitions that change 
the topological order by condensing bosonic vortex de- 
fects; the transitions we consider here were showrP^ to 
be dual (in a certain limit) to the TFIM. 

We begin by describing KZ scaling for these Ising tran- 
sitions, showing that the dynamics and scaling properties 
are exactly equivalent to those of the Z2 gauge theory, 
and identifying the coarsening degrees of freedom of the 
string net that is condensed in the starting topological 
phase. However, the mapping is not an equivalence. The 
full spectra of the models are quite different, which is 
reflected in the energetic splitting of different ground- 
state sectors. We then consider perturbations away from 
this limit, finding that as for the Z2 gauge theory, in 
the scaling limit these other perturbations do not alter 
the coarsening dynamics, but can be either irrelevant or 
dangerously irrelevant perturbations. 



A. Levin- Wen Hamiltonians with Ising transitions 

The SU(2)fc models we study live in a Hilbert space 
built from tensoring a finite set of spin variables on each 
link of a honeycomb lattice, 01 £ {0, |, 1, |}. These 
are analogous to the set of possible electric fluxes (of = 
±1) in the Z2 gauge theory. The idealized Levin- Wen 
Hamiltonians are similar in spirit to the toric code, and 
are constructed from a set of commuting projectors 



-KY,Vp-T m Y,Vs 



(22) 



where P represents a plaquette, and s a site. These mod- 
els can be viewed as deformations of lattice gauge theories 
with a continuous gauge group. They are 'deformed' in 
the sense that their representation theory is truncated, 
even though the gauge group is not discrete. In our con- 
text, a lattice 5/7(2) theory would have electric fluxes 
corresponding to all allowed spin values 0, 1/2, 1, while 
the models in question have a maximum spin ft/2. (We 
refer to the link spins as 'electric flux' though, more 
accurately, they are the representations of the lattice 
gauge/quantum group). Instead of describing our analy- 
sis for general values of ft we will now specialize to ft = 2 
in the interests of pedagogical simplicity and return to 
comment on the generalization subsequently. 

We now discuss the detailed form of the Hamiltonian 



olations of angular momentum conservation, analogous 
to the Gauss's law constraint G s = 1 in the Z2 theory. 
If the three links entering a vertex have spins i,j and /, 
angular momentum conservation requires I £ i x j. The 
rules for adding angular momentum have to be modified 
to be consistent with the truncation, however. For the 
SU(2) 2 model, the result ifpl 



(23) 



I (0,0,0) (0,1.1) 

otherwise 



(Hi) 



where it is understood that the eigenvalue of V s is in- 
dependent of interchanging the spins on the three links 
entering the vertex. 

The plaquette term Vp projects onto states in which 
P has no magnetic flux, and is written as a superposi- 
tion of "raising operators": Vp = ^(1 + J2t=i/2 a ^p)i 
where Bp raises all spins on the plaquette P by o in the 
truncated spin space. By "raising" , we mean a combi- 
nation of raising and lowering angular momenta in the 
truncated spin space. T> is the total quantum dimen- 
sion, equal to 2 here, while the coefficients a a depend on 
the quantum dimension^ of the spin representation o . In 
5/7(2)2, they are ao = 1, ai/ 2 = -v2, ai = 1. B P raises 
all spins in P by raising the spin on each link I £ dP. The 
action of Bf on a link I with spin i £ {0, 1/2, 1} is: 



Bj\0) = \i) 
^ 1/2 |l/2)(x|0)±|l) 
5/11/2) cx 1 1/2) 



B] /2 \l) cx |l/2) 
B}\1) cx |0) 



The numerical coefficients are chosen such that the am- 
plitude for creating any configuration with a 0-eigenvalue 
under V s is 0, ensuring that the vertex and plaquette pro- 
jectors commute. Their precise value is related to the 6j 
symbols of the quantum group 5/7(2)2, but we will not 
require their detailed form here. Interested readers can 
consult Ref. [20] for more details. 

As Vp and V s commute, the spectrum of the Hamil- 
tonian can be determined exactly. The ground state is a 
generalization (a string-net) of the loop gas ground state 
of the toric code, though there can be relative sign dif- 
ferences between terms in the Levin- Wen ground state 
wavefunction. As in the toric code, the excited eigen- 
states of (22) consist of "matter" excitations of energy 
Vmi and "vortex" excitations, of energy if. IntheSU(2) 2 
model there are anyonic spin-1/2 charges, fermionic spin- 
1 charges, and spin-1/2 or spin-1 vortices, both of which 
have bosonic statistics. The spectrum can be made to 
correspond exactly to that of the doubled SU(2)2 Chern- 
Simons theory. Accordingly the topological ground state 
degeneracy is knowrP^to be 9, as in the doubled Chern- 
Simons theory. 

We can drive a phase transition in our system by per- 
turbing the model ( 22 1 with transverse fields which create 



(22) for SU(2)2- The vertex projector V s penalizes vi- pairs of charges or vortices, as we did for the toric code by 
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adding erf and erf. The vortex excitations have bosonic 
statistics and hence transverse fields which create vortex 
pairs can drive a transition to a vortex condensed phase 
in which string-nets are confined. On the other hand, 
the analogue of the Higgs transition is not evident for 
our problem as both charges are non-bosonic. 

To drive the Ising transition that we are interested in, 
we add a transverse field which will condense spin-1 vor- 
tices. The Hamiltonian that we will tune through this 
transition is 



Hi 



SU(2), 



P s 

-kY: ^-rEl- 1 ^ 1 ( 24 ) 
p ^ z I 

where we have separated Vp into operators that "raise" 
spins by integer and half-integer amounts, and added a 
transverse field perturbation, T(— \) 2<Tl . The transverse 
field creates a pair of spin-1 vortices on the plaquettes 
adjacent to /, and has eigenvalue 1 on integer spin links, 
and —1 on half- integer spin links. Because the transverse 
field term squares to the identity (and all vortex creation 
operators commute), the vortices are Ising like. 

On every plaquette P and site s, the eigenvalues of 
Bp and V s are conserved, since these operators commute 
with -ffsu(2) 2 ■ Thus, we can consider the transition en- 
gendered by varying the ratio K/T in the subspace of the 
Hilbert space where the conditions 



I*) 



I*) 



(25) 



are always satisfied. In this subspace, the transition can 
be mapped onto the transition in the pure Z 2 gauge the- 
ory discussed in Sect. |II C[ and is therefore in the 3D 



Ising universality class as discussed in Ref. EH Here we 
give a different derivation of this result which focuses on 
the ground state wavefunctions and is better adapted to 
our purposes in this paper. 



To understand the mapping between H- 



SU(2). 



and Hy 



Q, notice first that the condition V s \ip) = \tp) ensures 
that we are always working in the "charge-free" sector 
where the (deformed) angular momentum is conserved 
at each vertex. This, together with Eq. (23), stipulates 



that we only need to consider configurations where the 
number of half-integer spins entering each vertex is even - 
or equivalently, configurations in which half-integer spins 
form closed loops. Similarly, in the absence of charge 
in the pure Z 2 theory §f§, the gauge constraint G s in 
Eq. (|]) ensures that links with electric flux (a x = — 1) 
form closed loops. The Levin- Wen transverse-field oper- 
ator (— l) 2ai assigns an energy penalty to the spin 1/2 
edges that form these loops, similar to the action of 
the transverse-field term erf on links with electric flux 
in the Z 2 theory. Thus both models describe a transition 
in which loops (of half-integer spin variables in the Levin- 
Wen case or o~ x — —1 variables in the Z 2 gauge theory 
case) become confined, and vortices become condensed as 
T/K increases. 



There is, however, a qualitative difference between the 

1/2 

operators Il;gap °f an d Bp i both of which change the 
number of loops in a given configuration. While erf sim- 

1/2 

ply flips the spin on the link Z, the operator B p maps a 
spin or 1 link to a spin 1/2 link, but a spin 1/2 link to 
a superposition of a link in the state and a link in the 
state 1. Thus one might worry that the two operators 
generate the same set of configurations (after identifying 
s = 0, 1 with a x = 1, and s — 1/2 with a x = —1), but 
with different statistical weights. 

We show in Appendix [B] that this is in fact not the 
case. Specifically, we prove that for any T, the ground 
state wave-function of either model can be expressed in 
the form 



I*) 



{1} 



{i}> 



(26) 



where {/} denotes a set of links on which cr; = 1/2 in 
the Levin- Wen model restricted to (25), or erf = —1 in 



the pure Z 2 gauge theory. Crucially, we find that /?{;} is 
the same for each set {1} in both models. Operators in 
the Levin- Wen model which commute with the conditions 



( 25 1 are either diagonal in the vortex basis, or diagonal in 
the spin basis and sensitive only to the spin on each edge 
modulo 1. The expectation value of any such operator is 
therefore identical to that of its Z 2 analogue, cementing 
the equivalence of the two models. 



We conclude that within the sub-sector ( 25 1 , the tran- 



sition is equivalent to that of the pure Z 2 gauge theory, 
and dual to that of the TFIM. It follows that our previ- 
ous discussion of string net coarsening, and the scaling 
of vortex creation operators, applies mutatis mutandis to 
the model at hand. 

Note, however, that the topological order of the initial 
and final phases of Eq. ( p4| is not the same as in the Z 2 
gauge theory; there are additional deconfined excitations 
on both sides of the transition. (In fact, the confined 
phase of -ffsu(2) 2 is a ^2 gauge theorySB). There must 
therefore be some operators in the Levin- Wen model 
whose behavior through the ramp is not captured by the 
mapping to the Ising gauge theory. 

To make this more explicit, we consider the fate of 
Wilson loop operators. In the Levin- Wen model there 
are two of these: 



W 1/2 (R,t-r) = (IJ^l (*) 



wee 



Wx{R,t;r) 



(27) 



\iec 



In order not to create vortices along the curve C, 
these raising operators must act with appropriate 
configuration-dependent complex coefficients, as dis- 
cussed in Refs. l20l and l52l Wi/ 2 {R, t; t) is clearly the 
analogue of the Wilson loop operator in the Z 2 gauge 
theory; its expectation value obeys a perimeter law in 
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the small T/K phase, and an area law in the large T/K 
phase, and its universal scaling in a linear ramp is given 
by Eq. (20). However, as the spin-1 variable remains 



deconfmed throughout the phase diagram, Wi(R, t; r) al- 
ways obeys a perimeter law. Its expectation remains con- 
stant in the scaling limit as the system passes through the 
critical point. 

A related effect can be seen in the amplification of the 
energy splitting between ground states in different topo- 
logical sectors discussed in the previous section. The 
SU(2)2 Levin-Wen model has 9 ground states on the 
torus. It can be showrP^that 6 of these obtain an energy 
cost of order L in the condensed phase, while one new 
ground state sector emerges, leaving a total of 4 topo- 
logically distinct ground states. These 4 correspond to 
the two possible eigenvalues of W" c (i?,i;T) along each 
non-contractible curve. 

Though we have primarily discussed the SU(2)2 Levin- 
Wen model, the main results apply to a large family of 
models in which there is an excitation that behaves like 
the Ising vortex^. Specifically, all the SU(2)fc models 
exhibit Ising transitions in which the half-integer spins 
(integer spins) can be mapped onto Z2 gauge configura- 
tions with a x — —1 (cr x = 1). They have two families of 
Wilson-line operators: the half-integral Wilson line oper- 
ators, which obey an area law in the confined phase, and 



scaling relations analogous to those of Eq. ( 20 ) ; and the 



integral Wilson line operators, whose expectation values 
do not depend on t, t and which remain perimeter law 
throughout the ramrJ^l. 

It is worth mentioning that the Ising transition we have 
discussed here is but one of a variety of confinin g tran- 
sitions that can be realized in Levin- Wen models ^ * 54 * 55 !. 
In the SU(2)2 model discussed above, for example, we 
could also add a transverse field term of the form cos irsi 
(which has eigenvalues (1,0,-1) for s/ = (0,1/2,1), re- 
spectively). This confines both spin- 1/2 and spin 1 la- 
bels, engendering a transition to a completely confined 
phase where both Wilson loop operators in Eq. ( 27 1 obey 



an area law. In this case the vortices that proliferate are 
not Ising-like, however, since the operator cos irsi does 
not square to 1. Very little is known about the critical 
theory in this case, and we expect that the transition is 
not in the 3D Ising universality class, so that the scaling 
functions and coarsening behavior will be fundamentally 
different from those of the Z2 gauge theory. 



B. Away from the pure Z2 limit 

Thus far, we did not have to concern ourselves with the 
other excitations in the SU(2)2 Levin Wen model as the 
their number was conserved in the ramp and we remain 
in the subspace (25) at all times. However, the other ex- 



citations, charges of spin 1 and 1/2, and vortices of spin 
1/2, will be created in a ramp if we perturb away from 
the limit of Eq. ( 24 ) by adding terms to the Hamiltonian 
which break the local conservation laws. In particular, 



the spin 1/2 charge has anyonic statistics relative to the 
condensing vortices, so that a finite density of these de- 
stroys confinement (these are analogous to the matter 
sources of the Ising gauge theory) . Once again, the KZ 
scaling limit saves the day. These spin 1/2 charges re- 
main gapped throughout the transition. Thus, as for the 
Ising gauge theory with matter, in the scaling limit we 
expect that the density of all of these excitations is van- 
ishingly small throughout the coarsening regime; terms 
violating the conservation of the spin- 1/2 charge at each 
vertex then act as dangerously irrelevant variables in the 
manner described in Section fVl 

Spin 1 charges and spin- 1/2 vortices, however, have 
bosonic statistics relative to the spin 1 vortex, and do not 
have analogues in the Z2 theory. Once again, for small 
perturbations which violate the exact local conservation 
of these excitations, they remain gapped throughout the 
transition and hence do not affect coarsening in the KZ 
scaling limit. Since a dilute density of such charges does 
not destroy confinement, however, we expect that they 
will not destroy coarsening even outside of the scaling 
limit. 



VIII. CONCLUDING REMARKS 

In this paper, we have initiated the study of the 
Kibble-Zurek problem for topologically ordered phases 
by studying the linear ramp across a transition that re- 
duces/breaks topological order and written down a scal- 
ing theory for it. Interestingly, unlike broken symmetry 
cases where it is natural to ramp from less to more order, 
here it is more natural to ramp from more to less order. 
The latter leads to our identification of the slow dynam- 
ics of string net coarsening much as the former leads to 
defect coarsening a la Kibble and Zurek. Of course, one 
can study the reverse protocol and the associated scaling 
although we have not done so here in the interests of not 
taxing the reader's patience unduly. 

The basic framework here can be easily generalized to 
other transitions out of topological phases; although for 
string-net coarsening to be visible, the gauge degrees of 
freedom must have a ready identification. Examples are 
transitions out of Z„ phases with n>3mci = 2 + l and 
with n>2ind = 3 + l. The Levin- Wen models also offer 
a "target rich" domain, although the analysis is likely to 
prove more complicated for more general condensation 
transitions. It will also be interesting to move to con- 
texts with conserved currents where one can study the 
temporal and spatial evolution of transport coefficients, 
such as the Hall conductance. 

Finally, for the statistical mechanically inclined, we 
would like to draw attention to our identification of 
gapped matter as a dangerously irrelevant variable in the 
dynamical KZ context. This is clearly a more general 
idea — e.g. irrelevant departures from intcgrability will 
be similarly dangerous — and it suggests that in the KZ 
problem, more couplings will be classified as such than 
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in the standard equilibrium analysis. 
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Recall that the growth law when £ is a function of time 



is (Eq. (17)) 



lco{t;r) ~ far) 



l/z d 



where in the last step, we have substituted the critical 
exponents of the (2+l)D TFIM, z = 1, z d = 2. Putting 
the pieces together, we see that t co <C ir =^ v < 1- 
This certainly holds at the 3D Ising critical point where 
1/^0.6. Thus, we conclude that coarsening described by 
Model C is indeed the correct long time asymptote for 
the KZ scaling functions in a linear ramp. 

Finally, we observe that all dynamical processes in the 
(2+l)D TFIM do not remain in equilibrium in the KZ 
ramp at late times. The scattering time between quasi- 
particles, t pv , grows as A 2 in this limit and is parametri- 
cally larger than ir- A hydrodynamical description, if it 
exists, is therefore more delicate than the case when the 
ramp is stopped at some t/t K = T . 



Appendix A: Scattering times in coarsening 

Here, we identify the dynamical process enabling coars- 
ening at late times alluded to in Sec. |III B| and justify 
that it remains in equilibrium during the KZ ramp. The 
criterion to remain in equilibrium is that the time-scale 
for such a process, t co , is parametrically smaller than the 
time-scale for the change in the transverse- field £r . Using 
(Tm — Tmc) = —t/T, we estimate t-p to be: 



M 



r 



Mc 



dFirf I dt 



= t. 



The system at late times has two kinds of excitations 
that are remnants of the paramagnetism at early times: 
1) The long domain walls of average size l co and 2) The 
bulk gapped quasi-particle excitations about each ferro- 
magnetically ordered state. As t/t K — > oo, the latter 
can be treated as classical particles. The average density 
of these particles and their momentum is essentially de- 
termined at t ~ t K and is fixed to be ~ 1/Z 2 and l/l K 
respectively. Their mass is determined by the gap A(i). 
The growing mass and the long inter-particle distances 
as compared to £ justify the classical particle approxi- 
mation. An average velocity of these particles can be 
determined as 



m 



l/l K 



where p is the average momentum and m the mass. The 
mechanism of coarsening proceeds through the transfer of 
energy between the long domain walls and these particles. 
To wit, the relevant time-scale t co is the scattering time 
between these particles and the wall: 



tc 



Ir 



Appendix B: Mapping of general SU(2)fc models to 
the Z2 gauge theory 

In this Appendix, we will discuss in more detail the 
mapping from the confining transition in the S f C/(2) 2 
Levin- Wen model to the pure Z 2 gauge theory |7]). As 
discussed earlier, the transition in question involves vary- 



ing K/T in Eq. (24), while restricting the Hilbert space 
to states with eigenvalue 1 under the vertex projector, 
and the integer part of the plaquette projector. In this 
subspace, links with half-integer spins form closed loops. 
The mapping to the pure gauge Z 2 theory involves map- 
ping half-integer (integer) spins to the presence (absence) 
ofZ 2 electric flux a x — —1 (+1). The transverse field op- 
erator (— l) 2a ' maps to af, and Bp to Bp = Jlieap °f ■ 
We will now show that the probability to be in any loop 
configuration is the same in both theories for every choice 
of K/T. 

We will begin with some notation. Let C denote the 
collection of links on the lattice that form closed loops, 
and (a!/ 2 (C)) the probability for a configuration in C. 
We will work here in the restricted Hilbert space of states 
for which 



B P \V) 



(Bl) 



and assume that our lattice has no boundary. To make 
the analogy to the Z 2 gauge theory, we also define 
the analogous operator, a x (C), whose expectation value 
gives the probability for a closed loop configuration of 
links with af = — 1. 

Our objective is to prove that, for every T and K, 
(a 1 / 2 (C)) = (a x (C )}. Since operators that commute 



with the conditions (Bl) arc cither diagonal in the spin 



basis and sensitive only to si mod 1, or diagonal in the 
vortex basis (dual to the basis of spin- 1/2 loops), this is 
sufficient to prove that their critical behavior is identical. 
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We will carry out the proof in two steps. First, we 
will show the equality for the two solvable points T = 
0, K > and K = 0,T > 0, where we can construct 
exactly the ground-states in both models. We will then 
use perturbation theory to argue that the result holds 
throughout the phase diagram. 



1. Equal weighting of loops in the ground states at 
the solvable points 

For K = 0, r > 0, the ground state has of = 1 in 
the Z 2 gauge theory, and oi £ {0, 1} for the Levin- Wen 
model. In this limit, for any C we have trivially that 
(ai/ 2 (C)) = (a x (C)) — and the result holds. 

Focusing on the opposite limit ( T = 0, K > 0), let us 
construct the exact ground states in the two models. We 
begin with the Z 2 gauge theory. Let |0) denote the state 
with o x = 1,t£ = 1 on all links and sites. This satisfies 
the Gauss law, but is not an eigenstate of the plaquette 
projector. To construct such an eigenstate, we take 



I* 



1 



TC) 



N P 

ee n 



B 



P u = 



n—1 



v 

N ^ 

{C} 



(B2) 

where *P„ runs over all possible distinct choices of n pla- 
quettes on the lattice, and N is a normalization. This 
sum generates all possible configurations C of loops with 
o x = —1, weighted equally (each configuration is in fact 
generated twice, since Yip Bp = 1). Since Bp = 1, 
b p z Ilp e *p n B P = Upe*P' n B P, where *P' n is *P„ with 
Pi either added (if it was not originally in the set) or 
deleted (if it was). It follows that 



1 



ee^ n b 



i 



Pe*i 



N P 

. EE II n >' 

n=l *P„ Pe*P„ 



and Y&tc) is a ground state. 

A similar construction can be used in the Levin- Wen 



models. Let Y$! e ) be a state satisfying (Bl|, with 07 an 
integer for every link I (careful inspection of these two 
conditions reveals that \^ e ) is a superposition of config- 
urations of closed spin-1 loops). Now consider: 



N P 



I* 



LW) 



^ee n 

v iv n =l *P„ P£*F n V 



(B3) 



We will show presently that 



Bi 



1/2V _ 



= 1 + B P Bp /2 B P = B P Bp /2 = Bp' 2 (B4) 



Using this fact, we have 



It follows that 



1/2 



(B5) 



and I^lw) is a ground state. 

Now, for any closed loop L, we can generate a configu- 
ration with 01 = 1/2 on all links in L, and no other links, 

1 /2 

by acting on \^ e ) with the product of B p on all pla- 
quettes inside the loop or all plaquettes outside the loop 
(these are the only such configurations on the right-hand 
side of Eq. (B3 )). In the Z 2 gauge theory the same holds 
for closed loops of o x ~ — 1. Hence given C, we have 

(a 1/2 (c)) = ( n -^ B p /2 ) i* e ) 

\pe*F(C) VZ J 

=^* e i n l(i+B P )\^) 



ps*p(C) 



1 

N 



(B6) 



where *P(C) contains either all plaquettes inside, or all 
plaquettes outside, the closed loops in configuration C, 



and the last equality is a result of imposing (Bl|. Thus 
for r = all possible configurations of spin-1/ 2 loops 
occur with equal probability in the ground state of the 
SU(2) 2 Levin- Wen model. 

We note that these results carry over directly to the 
more general case of an SU(2)fe Levin- Wen model, upon 
replacing spin- 1/2 (spin-1) with the set of all half- integer 
(integer) spins, and Bp 2 (Bp) with the sum of all half- 
integer (integer) spin-raising terms in the plaquette oper- 
ator (weighted by their respective quantum dimensions) . 



It remains to show that Eq. (B4| holds, which we will 
do for general k. We let Vp 2 — ^ J2a=i/2 3/2 a <jB P 



denote all half-integer raising terms in the plaquette op- 
erator, and Vp — ^ (l + J2 a =i 2 a<jB p ^J denote all 
integer terms. Following Levin and Wen, we choose the 
constant V such that Vp = \ (vp + Vp^ 2 ^j is a projec- 
tor. We then have 



(V P + Vp 12 ) 2 = (V F 



(V P 12 ) 2 



vi^vh + v\v 



1/2 



2(Vp + V 



1/2, 



Now, 
raise 

v P vp /2 

erators 



(V]/ 2 ) 2 and (Vp) 2 both contain only terms that 
the spins in P by an integer amount, while 



Vy 2 Vp contains only half-integral raising op- 



Pp. 



It follows that 



(Vp 12 ) 2 = 2Vp , 



VpVp 1 = Vp 



1/2 



(B7) 



B 



1/2 
p, 



n 



B 



1/2 



Pe* 




Pi € *P„ 
Pi i *P„ 



We also have 



(Vp + Vp 2 )(Vp - V 



1/2, 



(B8) 
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since it can be showrFU that (Vp + Vp 2 ) projects onto Thus is a state with vortices on each plaquette 

flux-free states, while (Vp — Vp^ 2 ) projects onto states 
with an Ising vortex. It follows that 



{V 1J2 )2 



(B9) 



This also implies that Vp is a projector, a nd t hus that 
the eigenvalues of Vp are and 1. (From Eq (B9 ) and the 
fact that Vp is a projector, it follows that the eigenval- 
ues of V X p 2 are 0,±1; when restricted to configurations 
where Vp\^) = \^f), they are ±1, as one expects from 

1/2 

the correspondence of V P to the plaquette term of the 
toric code.) 



with an odd number of edges in the set of links {/}. It 
also follows that choices of {/} which differ by a product 
n ;gC * hi, where C* is a set of closed curves on the dual 
lattice, create identical excited states, as rLec* ^zl^o) = 
|^o)- Finally, we have 



(Bll) 



In other words, the inner product is 1 if the combination 
of the two sets {1} and {/'} of links forms a set of closed 
curves on the dual lattice, so that ^{j}) and ^{p}) have 
vortices on the same plaquettes. Similarly, we may com- 
pute matrix elements of the Hamiltonian within these 
excited states via: 



2. Away from the solvable points 

Next, we wish to show that the result of the previous 
section holds true throughout the phase diagram. One 
way to do this is to invoke the result of Ref. 48, where 
it was shown that within the subspace of states satisfy- 
ing (Bl I, the SU(2)fc Levin- Wen models are exactly dual 
to the transverse-field Ising model. We can identify all 
states in this Hilbert space by the configuration of dual 
Ising spins (together with their topological ground-state 
sector, in periodic boundary conditions). The duality 
relation — which also holds for the Z 2 gauge theory - 
ensures that the probability amplitude to be in a given 
vortex configuration is identical in both models. The 
physical operators in this Ising subspace are either diag- 
onal in the vortex (or dual Ising spin) basis, or diagonal 
in the basis of spin-1/2 loops. (These are precisely the 
operators that do not cause violations of (Bl), and can- 
not distinguish between edges of spin and spin 1). It 
follows that all expectation values of such operators - 
including (a 1 / 2 (C)} — must also be identical to their Z 2 
analogues (such as (a x (C))). 

Here we will take an alternative, perturbative approach 
to prove the desired result. We will begin at an arbitrary 
point in the deconfined phase, and consider construct- 
ing the wave-functions in both theories to some finite 
order in perturbation theory. These wave functions are 
linear combinations of the unperturbed (Levin- Wen or 
toric code) ground state, together with excited states of 
the form 



1*0}) 



n 

Z6{Z} 



(BIO) 



where we have defined the transverse field operator hi = 
af for the toric code, and (-T) 2s ' for the Levin-Wen 
model, and \^o) denotes the unperturbed ground state. 
If l\ and l 2 are two links bordering plaquette P, we have 



.B f 



-B P at 



J h u h> 



p ~ 
Bp] = 



-B l J 2 (-iy°^^ 



\2s. 



(-l) a «i(-l) 2 



, B 



1/2 



= 



<*«i n * 

ieO"} 



^{/}u{/'}u{i"},C* 



(B12) 



The crucial point is that for any choice of {I}, {I'}, {I"}, 
these matrix elements are identical in both models. Since 
the weight of each unperturbed excited state in the ex- 
act ground state can be constructed perturbatively using 
only matrix elements of this form, it follows that 



i*> = E^}i*w> 



(B13) 



{'} 



with /3{;j the same for each set {1} in both models. 

Finally, we observe that ai/ 2 (C) and a x (C) both have 
the form 



(B14) 



and, in particular, commute with hi on every link. (Here 
a v = Cii /2(C) , a x (C) as appropriate). This, together 
with the relation ( |B12[ ), implies that 



01,0'} 



n h ' a » n &ii*°> 

ie{i} ie{i'} 



in both models. We have already shown that the co- 
efficients /3m are the same, and the possible choices of 
{I}, {I'} on which the 5 function has support are a geo- 
metric property of the lattice. Invoking the result of the 
previous subsection, we can thus conclude that for all T 
in the deconfined phase, 



(ai/ 2 (C)) r - K(C))j 



(B15) 



Our derivation has implicitly relied on the fact that 
we can construct the exact ground state perturbatively, 
starting from the ground state of the toric code or Levin- 
Wen solvable point. Thus the above argument fails at 
the critical point, and in the phase where T/K is large. 
In this regime, however, we may make essentially the 
same argument, by replacing hi with the plaquette opera- 
tor, and I^lw), |*tc) (denoted by |* ) m the derivation 
above) with |^ e ) and |0) respectively. In this case, the 
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basis of excited states generated will be an eigenstate of 
a x (toric code) or (~l) 2s (Levin-Wen). There is no need 
to define an analogue of C* , since if two distinct products 
of plaquette projectors produce the same loop configura- 
tion state, then their product is the identity operator. 



In each phase, we can thus argue that Eq. (B15 ) holds 



to arbitrary order in perturbation theory. It follows that 
as the phase transition is second order, it must also hold 
at the critical point, proving the result. 
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